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I. INTRODUCTION 

In the vicinity of the glass transition point, the dynamics of supercooled liquids becomes extremely slow [1-3]. 
The dynamics of the glass transition attracted much attention over the years. Among many theoretical approaches, 
the mode-coupling theory (MCT) is one of the most successful ones which can be "derived" from the first principle 
equation, and explains many aspects of observations in experiments and simulations, such as multi-step relaxation 
processes and the Debyc- Waller parameter [4-7]. 

In spite of such advantages of MCT, there are some controversial points for the validity of the standard MCT 
(SMCT). Indeed, SMCT predicts the existence of the ergodic-nonergodic (ENE) transition, where the system becomes 
nonergodic below a critical temperature or above a critical density, while real systems are still ergodic in experiments 
and simulations at low temperature or high density. Furthermore, SMCT predicts an algebraic divergence of the 
viscosity at the critical point of ENE transition, but the viscosity for real supercooled liquids obeys the Vogel-Fulcher 
law near the glass transition point and the Vogel-Fulcher temperature is lower than the critical temperature of ENE 
transition. To overcome these difficulties of SMCT, many investigations have been carried out [8-29]. The failures 
of SMCT may be originated from the decoupling approximation of a four-point correlation function. In fact, Mayer 
et al. [21], introduced a toy model which does not have any spatial degree of freedom, and demonstrated that the 
ergodicity of the system at the low temperature is recovered when they include higher-order correlations, while there 
exists ENE transition within the framework of the decoupling approximation. It suggests that we should not adopt 
the decoupling approximation, but use an approximation which contains higher-order correlations. However, the 
systematic improvement of the approximation is difficult within the conventional framework by using the projection 
operator technique. 

The field theoretical approach is the promising method which can systematically improve the approximations. There 
is another advantage of the field theory in which we can discuss the response function and the fluctuation-dissipation 
relations (FDR). Following Martin-Siggia-Rose (MSR) method [30], we can construct an action by the introduction 
of conjugate fields, for a set of nonlinear Langevin equations, and can use the perturbative expansion. The current 
situation for this approach, however, is confusing. Indeed, among many field theoretical investigations [8, 10, 12, 20, 
23, 27-29], only a few papers have succeeded to derive SMCT in the lowest order perturbation from the nonlinear 
Langevin equations. One of main difficulties lies in the violation of FDR in each order of naive perturbative expansions 
of the set of nonlinear Langevin equations, as indicated by Miyazaki and Reichman [20]. 

In order to recover FDR-preservation at each order of the perturbation, recently, Andreanov, Biroli and Lefevre 
(ABL) [23] indicated the importance of time reversal symmetry of the action, and introduced some additional field 
variables. Indeed, ABL demonstrated that we can construct a FDR-preserving field theory, starting from the nonlinear 
Langevin equations which contain both the Dean-Kawasaki equation and the fluctuating nonlinear hydrodynamic 
(FNH) equations. Kim and Kawasaki [28] further improved ABL method and they derived a mode coupling equation, 
similar to SMCT, from Dean-Kawasaki equation [11, 31] in the first-loop order via the irreducible memory functional 
approach which to be essential for treating the dynamics of the dissipative systems such as the interacting Brownian 
particle system. 
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On the other hand, Das and Mazenko [8] pubHshed a pioneer paper on the field theoretic approach of FNH. They 
suggested the existence of the cutoff mechanism in which the system is always ergodic even at low temperature. 
Later, Schmitz, Dufty and De (SDD) [10] reached the same conclusion as that by Das and Mazenko from a concise 
discussion, though they destroyed the Galilean invariance of FNH equations. On the other hand, Kawasaki [11] 
suggested that FNH equations reduce to Dean-Kawasaki equation in the long time limit. Furthermore, ABL [23] 
suggested the existence of ENE transition in FNH, and indicated that the calculation by Das and Mazenko breaks 
FDR-preservation. Moreover, Gates and Ramaswamy [22] indicated that the calculation by Das and Mazenko violates 
the momentum preservation. Das and Mazenko [32], however, responded that the indications by ABL and by Gates 
and Ramaswamy are not the fatal errors of Das and Mazenko [8], but contain some misleading arguments. Thus, we 
are still in a confusing situation for the application of the field theory to the glass transition, and cannot conclude 
whether ENE transition exists in FNH. 

In this paper, we apply the method developed by Kim and Kawasaki [28] to FNH to clarify the current situation 
of the FDR-preserving field theoretical approach to the glass transition. The organization of this paper is as follows. 
In the next section, we introduce FNH which describes the time evolutions of the density field and the momentum 
field, agitated by the fluctuating random force, for compressible fluids. This set of equations is equivalent to that 
used by Das and Mazenko [8] and ABL [23]. In the former half of Section HI, we make an action invariant under 
the time-reversal transformation. In order to keep the linearity of the time-reversal transformation, we introduce 
some additional variables and their conjugate fields. This linearity of the time-reversal transformation makes FDR- 
preserving field theory possible. We also introduce a complete set of Schwinger-Dyson equations of our problem, and 
summarize some identities used for the perturbative calculation in the latter half of Section III. Section IV is the main 
part of our paper, in which we explain the detailed calculations of perturbative expansion within the first-loop order, 
under the assumption that the correlations including momentum can be ignored in the long time limit. Within this 
approximation, we predict the existence of ENE transition, and reach an equivalent equation obtained from SMGT. 
In the last section, we discuss the validity of our assumptions used in this paper, and compare our results with 
others. We also summarize our results. In Appendix A, we introduce the details of the time-reversal transformation 
and some relevant relations derived from the time- reversal transformation. In Appendix B, we present the details of 
the calculation for one component of the Schwinger-Dyson equation. In Appendix G, we show some relations for the 
equal-time correlations and the self-energies. In Appendix D, we write down the explicit expressions for all three-point 
vertex functions. 



II. FLUCTUATING NONLINEAR HYDRODYNAMICS 



In this section, we briefly summarize our basic equations, FNH, and MSR action [30]. The argument in this section 
is parallel to those in the previous studies [8, 23]. 

Let us describe a system of supercooled liquids in terms of a set of equations for the density field p{r, t) and the 
momentum field g{r^ t). For the continuity equation of momentum, we employ Navicr-Stokcs equation for compressible 
fiuids supplemented with the osmotic pressure induced by the density fluctuation, and the noise caused by the fast 
fluctuations. In order to keep the analysis simple, we ignore the fluctuations of energy [33] as assumed by Das and 
Mazenko [8], SDD [10] and ABL [23]. 

The time evolutions of the collective variables p and g, which we call FNH equations, are given by [8, 23] 

dtp = -V-g, (1) 

Ot9a = -p\'a-^ Lai3 h T^a . (2) 

dp ^ p p 

Here, rja is the Gaussian white noise with zero mean, which satisfies 

(77„ (r, t) T^p (r', t')) = 2TL^p5 (r ~ r') 5 {t - t') , (3) 
where T is temperature and L^p is the operator tensor acting on any field variables h (r) as 

Lc^ph (r) = - Qv„V/3 + S^pV"^^ + CoVaV^jj h (r) , (4) 

with the shear viscosity /iq and the bulk viscosity Co- We note that here and after the Boltzmann's constant set to 
unity. In this paper, the Greek indices, such as a, are used for the spatial components, and Einstein's rule such as 
gaQa = Yl^a=i 9a ^'1^0 adopted. The effective free-energy functional F = Fk + Fjj consists of the kinetic part Fk 



3 



and the potential part Fu as 



Fk =1 I dr 
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(5) 
(6) 



where m is the mass of a particle and c (r) is the direct correlation function [34]. The potential part Fjj of the effective 
free-energy functional has the same form with Ramakrishnan-Yussouff form [35]. Here, dp (r, t) = p (r, t) — po is the 
local density fluctuation around the mean density pq. From the relations (5) and (6), we can rewrite (1) and (2) as 



dtp = -V p 
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where we have used p\/a{SFK/Sp) = —g/jVaiSF/Sgp). 

In general, it is impossible to solve the set of nonlinear partial differential equations (l)-(6). In this paper, wc adopt 
MSR field theory [30]. Let us derive the MSR action. Because the collective variables p and g satisfy the dynamic 
equations (7) and (8), the average of an observable A[p,g] is expressed as 



(A) 



I Dp'Dg'A[p\g']6{p'-p)5{g'-g: 
= J DpDg J{p, g)A[p, g] / 6 (^dtp + V (^p 
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where J{p,g) is the Jacobian. As written in Ref. [36], the Jacobian J{p,g) can be independent of both p and g when 
we employ the Ito interpretation. When we replace the delta functions by the functional integrals of the conjugate 
fields p and ga, the average of A in eq. (9) can be rewritten as 

^ ^ DpDgDpDgA[p,g] 
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where Zq is the normalization constant. By means of eq. (3), the average of A is given by 
(A) ^ 1-J DpDgDpDgA[p,g] 

drdt 



X exp 



X (^cxp 



SF 



SF 



dtga + P^a-r. h ^[) ga-r. 
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= DpDgDpDgA[p,g]e^iP'P'a'^\ 
where the MSR action S[p,p,g,g] is defined by 
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S[p,p,g,g] = J drdt -pj 



SF 
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III. THE CONSTRUCTION OF A FDR-PRESERVING FIELD THEORY 
A. The time-reversal symmetry in the action 

In order to construct a FDR-preserving field theory, it is necessary to introduce some new variables, and the linear 
time reversal transformation, which makes the MSR action invariant. It is easy to check that the action (12) is 
invariant under the time reversal transformation [23] 

1 SF 

t^-t, p^p, p^-p+--z-, 
Here, we adopt the method developed by Kim and Kawasaki [28] in which they introduced the new variable 



Sp 



A*5p, (14) 



where A*5p represents the linear part oi 5Fu / 5p on 5 p. Kim and Kawasaki [28] confirmed that the density correlation 
function of Dean-Kawasaki equation for the non-interacting case satisfies the diffusion equation nonperturbatively. 
Furthermore, they concluded that the nonergodic parameter is same as that of SMCT. 

Following the idea of Kim and Kawasaki [28], to eliminate the nonlinearity of time reversal transformation of (13), 
we introduce the new variables 9 and u 

e=^^-^-K*5p, (15) 
I SF I 

Va^--. 7T^9a-, (16) 

J oga J PO 
where the operator K acts on any function h as 

K*h{r) = — [ dr' Is (r - r') - ^c(r - r')] h(r') . (17) 
mpo J I mi 

It should be noted that the right-hand side (RHS) of eqs. (15) and (16) do not include the zeroth and the first order 
of Sp and g. The choices of eqs. (15) and (16) differ from those by ABL [23]. The implication of the difference will 
be discussed in Section V. 

As the result of the introduction of 9 and j/q, the action (12) can be rewritten as 

S[ip]= J drdt -p{dtp + Va{p{p^^ga+Tv^))}-g^\dt9a + TpVc.{K*5p + e)+L^fi{p^^gfi + Tvp) 



+ V/3(5q(Po ^5/3 + Tvp)) + gii^aiPo ^9P + Tl^p) \ + TgaLapgii ~ 9{e - fg) - Va {Va - fuj 



(18) 



where we have introduced 



1 SF 

fg{Sp,g) = -—-K^Sp, (19) 

UASp,g) ^ ^^--!-9a. (20) 
^ Sga J PQ 

We have also used the abbreviation of a set of the field variables xp'^ = (^Sp, p, 9, 9, g, g, v, i>j . Here, the time reversal 
transformation which makes the action (18) invariant, is given by 

t ~t, p p, p ^ -p + 9 + K * Sp, 

ga -~9a, ga ,9a — ^^a ^ 5q: 9^9, 

Tpo 

9^9 + dtP, Va^~Va, Va ~Va - dtga- (21) 
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We, thus, can construct a FDR-prcscrving field theory, due to the hnearity of the time reversal transformation (21). 
As in the usual cases, let us split the action (18) into the Gaussian part Sg which represents bilinear terms of the field 
variables and the non-Gaussian part Sng as 

S[xP]^Sg[xP]+Sngm, (22) 

where 

Sg[xp] = j drdt\^-p {dtp + V«ff« + TpoV^Ug ] - 9o,\dtga + Tp^V^K * Sp + TpoV„6> 



(23) 



and 



Sng [tp] = J drdt I -p | v„ {Sp {po^gg^ + Ti^o,)) } - .9„ {T(5pV„ {K *5p + 6) 

+ "^liidciPo^gp + Tvp))+gpV^{p^^gp + Tz.^})} + efg{6p,g) + (^p,g)|. (24) 

Note that we present some relations in time-reversal symmetry of this action in Appendix A. 
It should be noted that the continuity equation (7) can be rewritten as 

dfP = -VaipiTVa + Pa^ga)) 

= ■g-TpnVa-l'a-^a{5p{TVa+ Po^ga))i (25) 

where we have used cq. (16). From eqs. (1) and (25) we obtain the identity 

TpciVaVa + Va{5p{TVa + Po^ga)) = 0. (26) 

Therefore, the sum of the underlined terms in eqs. (23) and (24) should be zero. However, each the underlined term 
is included in the Gaussian part (23) or the non-Gaussian part (24) separately. To satisfy the action invariant under 
the time-reversal transformation in 'each' part, we should keep these terms in the calculation. 

B. The exact results of the Schwinger-Dyson equation 

In this subsection, we derive a set of closed equations of two-point correlation function. Let us express the two-point 
correlation function in the matrix form as 

G{r-r',t-t')^{^{r,t)i,^{r',t')), (27) 

and its ipTp' component is represented by 

G^^, {r -r',t- t') = (V (r, t) i,' (r', t')), (28) 

where -0 or V'' is the one of the components of We note that here and after we adopt the simple notation ip ~ p 
to represent the contribution from 6p. With the aid of the Fourier transform of h (r) 

/rik 
-—^e^'^-hik), (29) 

and the action (22), we obtain the Schwinger-Dyson (SD) equation 

[Go"' • G] (fc, t) - p • G] (fc, t)^\S{t), (30) 

where Z and I are, respectively, the self-energy matrix and the unit matrix. Here, the free propagator matrix Go 
satisfies 

SgW dXidX2V{Xi)Q.a-\X^-X2)ip{X2), (31) 
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where we have used the abbreviations as Xi = (r^, ti) with i = 1,2. Wc note that the SD equation (30) is an equation 
for 16 X 16 matrices. 

Here, we indicate that Gpp{k, 0) is related to the static structure factor S (k) as 

S{k) = —Gpp{k,0). (32) 
mpo 

From the definition of the direct correlation function [34], eq. (32) can be rewritten as 

Sik)=(l-^c{k)Y' = —K-\k). (33) 
V m / mpo 

Let us explicitly write some components of the SD equation 

9tGp0 (fc, t) + ika,Gg^4, (fc, t) + iTpokaG^^^, (fc, t) = Fp^ (fc, t) , (34) 

Ge4, (fc, t) + ^gg (fe, 0) Gp^ (fc, t) - Fg^ (fc, t) , (35) 

(fc,O)Gg,0(fe,i) -^^^„0(fe,i), (37) 

where is the one of the components of (fp" = (5p, 0, g, i/) . Here, eqs. (34), (35), (36) and (37) are obtained from pcf), 
04), ga4> and Dacj) components of the SD equation, respectively. The derivation of eq. (35) is shown in APPENDIX B 
as one example. The derivation of the other equations is parallel to that for eq. (35). We express the RHS of eqs. (34)- 

(37) as a unified form F^^ where 4> is the one of the components of cf)^ = {j>, §,g, i>j . From the parallel argument to 
derive Fg^, i.e. RHS of (B5), F^^ satisfies 

P^^f^^t) = I rfs|-E^-(fe,t-s)(X(fe)G,^(fc,5) + Ge0(fc,s)) + S^^(fc,t-s)9sGp0(fc,5) 



{k,t-.s)d,Gg^^{k,s)\. (38) 



These components of SD equation are so complicated because of the self-energies. However, we can simplify the 
components of the SD equation with the aid of some exact relations. First, we note that the equal-time correlation 
functions satisfy (see (C4) and (C8) in Appendix C) 

EQg(fe,0) = 0, and (fe, 0) 0. (39) 

With the aid of eqs. (39) , eqs. (35) and (37) are, respectively, simplified as 

Ge^{k,t)^Fg^{k,t), (40) 

G,„0(fe,i) = Fe„0(fc,t). (41) 
Second, from eqs. (23) and (24), there is the following identity 

,6p{r,t)' 

This identity can be expressed by the explicit form 

{[dtp + y-g + TpoV ■ V ■ {5p{Tu + p^^g))] (r, i)0(r', t')) = 0. (43) 
With the aide of the identity (26), the Fourier transform of this equation becomes 

9tGp0(fc,<) + zA;„Gg„0(fe,O = O. (44) 

This equation implies that our SD equation preserves the mass conservation law. We also derive the identity by the 
substitution of (44) from (34) 

Fp^ (fe, t) = ik^TpoFo^^ (fe, t) . (45) 



r',t'))=0. (42) 
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IV. PERTURBATION IN THE FIRST-LOOP ORDER 

111 this section, we develop the perturbative calculation of the SD equation within the first-loop order approximation. 
When we assume that the correlations include the momentum decay fast enough, we can obtain an equation for the 
non-ergodic parameter in the long time limit. 

From eqs. (36), (40), (41) and (44). we, thus, obtain the time evolution of the density correlation function as 

d^Gpp (fe, t) + p^'LdtGpp (fc, t) + Tpok^K{k)Gpp (fc, t) 
= -Tpnk^Fg^ (fc, t) - ik^Fg^p (fc, t) + iTLk^Fp^p {k, t) , (46) 

where L = SapLap. This is a remarkable result that the left-hand side (LHS) of eq. (46) is equivalent to SMCT 
without memory functions when we omit the terms which include the self-energies. However, this equation is quite 
complicated, because the self-energies are included in Fg^, Fg^p and Fp^p. Therefore, we restrict our interest to the 
calculation of the self-energies in the first-loop order perturbation in the latter part of this paper. 
In the first-loop order perturbation, the self-energy S^^^, is expressed as 

I]^^^,(Xi,XO = i J2 J dX^dX^dX^dX;, 

^0102 03 (^1' ^2, ^3) ^4,[4,'^4,'^ {^11^2^ ^3) ^02 0^ {X2 ~ X2) G^g^!^ (X3 - X3) , (47) 

where or (jy'i is the one of the components of (/>, and 0i or 0^ is the one of the components of (f>. Here the three-point 
vertex V? , , is defined by 



F;,,,,3 ix^,x,,x,) = scpux^)lMi%<l>3{x,y 



We list the all three-point vertices which include cp in APPENDIX D. Note that there are no four-point correlation 
functions including both of (f)i and (f)[. 

Within the first-loop order approximation, F^^ in eq. (38) is reduced to 

C^'*) - I dsl^j:^g{k,t~s)dsGpAk,s)-iTpo)-'^^-Jk,t~s)Gg^^ik,s) 

-E^g^ (fe, t - s) Fo,^ {k, s) + E^.^ (fc, t - s) {dsGg^^ik, s) + ik^Tpo{K{k)Gp^{k, s) + Gg^{k, s))}| 

- J^ds |e^^ (fc, t - s) dsGp^ (k, s) - (rpo)-'S^g„ (fc, t - s) Gg^4, {k, s) 



'^4.i>Jf^^i- s)Po^^c,f}Gg^^{k,s)y 



(49) 



We have used (36), (40) and (41) and eliminate higher order contributions to obtain the final expression. We have 
also used the (41) and the relation 

E^^(fc,t)~-zfc„rpo%„(fe,t), (50) 

which is only valid in the first-loop order, for the first equality in eq. (49). Let us show the expression (50). The 
self-energy E^^ should contain the vertex Vppg^ or Vppi,^ . Among these two vertices, the vertex Vppi,^ is irrelevant 
within the first-loop order calculation. Indeed, the self-energy with the vertex Vpp^^ should contain the propagator 
G^^0 which is equal to F^^^f, from eq. (41) and is the first loop order. Thus, eventually, the self-energy becomes higher 
order correction as in Fig. 1. On the other hand, from eqs. (Dl) and (D9), the vertices Vppg^^ and V£>^pg^ satisfy the 
relation 

Vppg^, (Xi, X2, X-s) = Tpo\7r,c.V,^pg, (Xi, X2, X3). (51) 

Thus, we reach the relation (50). 

Let us calculate some typical terms, such as Fg^ (fc, t), which appear on the RHS of eq. (46) in the long time behavior 
under the first-loop order approximation. For simplicity, we assume that the correlations include the momentum decay 
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FIG. 1: One of the diagrams of the self-energy E^^ which is produced by the vertex Vpp^^. From the eq. (41), the first-loop 
order self-energy in the left-hand side of the figure can be treated as second-loop order as figured in the right-hand side of the 
figure. 



fast enough to be negligible in the long time behavior. For this purpose, at first, we calculate Egg (fc, t). Among the 
three-point vertex functions listed in APPENDIX D, there are only two vertices (D3) and (D4), which include 9. 
Substituting (D3) and (D4) into (47) with ~ 4>'i ~ 0, the expression of Egg (fc, t) at the first-loop order is given by 



+ TFZr^ (<^9»P (9' ^) <^ScP (fc -Q,t) + Gpg^ (q, t) Gpg^ (fc - q, t)) 



1 

Tmp^ 
2 



-j;2^G,„,,(9,t)G,„g,(fc-q,t)j. (52) 

In the limit t oo, thus, the first term of Fg^ (fc, t) in cq. (49) can be approximated by 

/ dsEgg{k,t- s)dsGpp{k,s) ~ Sgg(fc,t) / dsdsGpp{k,s) 
Jo Jo 

Here, the last expression is obtained from the assumption that the correlations include momentum decay fast enough. 
Similarly, with the aid of (47) and (D3)-(D8), Eg. {k,t) at the first-loop order calculation reduces to 

^eg^ C^' - - [ 77^ (k {q) Gpp (q, t) + Gap (g, t)) Gpp (fc - g, t) 

J (27r) 'tiPq \ / 

= I 7^<i4k (l)Gpp{Q,t) + Fg {q,t))Gpp{k - q,t) 

' f 'i"^' ^) ^PP ~ 9' 

J ( 27r ) 

kpqaK (q) Gpp {q, t) Gpp (fc - q,t) , (54) 



wPo J (27r) 
iTkfj f dq 



mk^pl J (27r)' 

in the limit t — > oo. The first equality in (54) comes from the assumption that the correlations include momentum 
decay fast enough. For the second equality in (54) we have used eq. (40). To obtain the third equality in (54) we have 
ignored the contribution from Fg^. This simplification can be justified at the first-loop order approximation, because 
Fg^ is the first or the above loop order function. To obtain the last expression in (54), we have used the fact that the 
density correlation function depends on time and the absolute value of the wave vector. From cq. (54), the second 
term of Fg^ {k,t) in cq. (49) becomes 

1 /•* 

- L 3 / T^^kaqaK {q)Gpp{q,t)Gpp{k- q,t) [ dskpGg^p(k,s) 
fni^ Po J (27r) Jo 

= , 2 3 / ^^,3 kaQaK (q) Gpp (q, t) Gpp (fc - q, t) ( dsdsGpp{k,s) 
™fc Po J (27r) Jo 

= ^ Gpp{k,t)-Gpp{k,Q) I (55) 

Po J (27r) 
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where we have used eq. (44) in the second cquaUty. The last term of F^^ in eq. (49) is zero because (A;,t) and 
Ggi^p{k,t) are zero in the long time limit. 

Thus, we obtain the expression for Fx {k,t) in eq. (49) from the eqs. (53) and (55) at the first-loop order as 



^o,., o ... .„t)(Gpp(fc,t)-Gpp(fc,0)), (56) 

(27r) 

in the limit t oo. Similarly, we evaluate ikaFg^p {k,t) and F^^p {k,t) within the first-loop order as 

ikaFg^p{k,t) = I I ^kaqaK {q) + {k^qaK {q) + ka{ka- qa)K {k - q)f\ 

J (27r) 1^ mpf^ 2fc^/9o J 

xGpp (q, t) Gpp (k - q, t) [Gpp (fc, t) - Gpp (fc, 0)) , (57) 

and 

Fo^p{k,t) = 0, (58) 
in the limit < — > cxo. 

From the above expressions for Fg^ {k,t), Fg^p {k,t) and Fp^p {k,t), we can evaluate RHS of eq. (46) in the limit 
t — > cx). We note that the first and the second terms on the LHS of eq. (46) are zero in the long time limit, because 
the time derivatives of the density correlation functions should be zero in such a region. Therefore, the self-consistent 
equation of the noncrgodic parameter /(fe), which is defined by 

/(fc)^ hm^^^^iiM, (59) 
t^oo Gpp (fc,0) 

is obtained as 

where 

with 

Vk^q = ka{qaC{q)+PaC{p)), (62) 

where we have used p = k — q and the static structure factor (32). This set of self-consistent equations (59)-(62) for 
the nonergodic parameter is equivalent to that in SMCT. 



V. DISCUSSION AND CONCLUSION 



A. Discussion 



In this paper, we formulate the FDR-preserving field theory for FNH. The obtaining SMCT under the first-loop 
order approximation in the previous section, is a successful first step to construct a correct theory beyond SMCT. 
However, we still have some unclear points in the analysis. Let us discuss such unclear points through the comparison 
with other field theoretical approaches. 

To analyze the SD equation, we employ the important assumption that the correlations include the momentum are 
negligibly small in the long time region. This assumption is crucial to discuss whether ENE transition exists. Indeed, 
if we assume that the contributions from momentum are negligible, FNH equations are reduced to the Dean-Kawasaki 
equation, as indicated by Kawasaki [11]. On the other hand, there are several indirect evidences for the justification of 
this approximation. First, we note that a numerical simulation exhibits fast relaxations of the correlations including 
the momentum [37]. Second, it is known that the density-density correlation Gpp is connected only to the correlation 
in the longitudinal mode k^Gg^^, which is proportional to dtGp^. Since the system is almost stationary, any terms 
including the time derivative are small. Thus, the contributions from the correlation including the momentum can be 
ignored in the slow dynamics in the motion of the density field. 
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paper 


(i) model 


(ii) FDR 


(iii) approximations 


(iv)iionergodic 
parameter 


DM [8] 


FNH 


non- preserve 


non-perturbative 





SDD [10] 


FNH witli violation of 
Galilean invariance 


preserve 


first-loop 

Gg^ 


order, 





ABL [23] 


Dean-Kawasaki and 
FNH 


preserve 


first-loop 

Gg4, ^ 


order. 


1 


KK [28] 


Dean-Kawasaki 


preserve 


first-loop order 


fsMCrik) 


This work 


FNH 


preserve 


first-loop 

Ggtf, —> 


order. 


fsMCrik) 



TABLE I: The comparison of our results with Das and Mazenko (DM) [8], SDD [10] and ABL [23] and Kim and Kawasaki 
(KK) [28]. We list the results on the four points: (i) the used model, (ii) whether FDR is preserved in the perturbation, (iii) 
the used approximation to derive the density correlation function in the limit t ^ oo, and (iv) the behavior of the nonergodic 
parameter f{k). Here, the expression Gg^ —f means that the correlations include the momentum becomes zero in the long 
time limit. fsMCT{k) means the nonergodic parameter is equivalent to that of SMCT. 



Let us compare our results with those obtained in other field theoretical researches, taking into account following 
four aspects, (i) Is the basic model adequate? (ii) Is the analysis FDR preserved? (iii) Are the approximations to 
analyze the density correlation function valid in the limit t —^ oo? (iv) What is the behavior of the density correlation 
function in the limit t ^ oo? The results of the comparison are summarized in Table I. 

First; we compare ours with ABL [23]. As can be seen in Table I, they predicted that the nonergodic parameter is 
unity, which is independent of the wave number. This result is clearly in contrast to the observations in experiments 
and simulations. This result might be caused by the definition of the new set of additional fields. Indeed, we have 
introduced 9 in eq. (15) and v in eq. (16), while ABL [23] 

SF dF 

Oabl = and v^^abl = — • (63) 

dp dga 

These new variables include the linear terms 5p and g, but our additional variables 9 and v do no include the 
linear terms. As a result, the order of the correlations, which include the new variables 9 and f, differ from ours. 
Indeed, Gg^ and G^^^ are the first or the above loop order in eq. (40) and eq. (41), while Ge^ and G^^^ include the 
tree diagrams in the calculation of ABL. Therefore, we suggest that eq. (63) is not appropriate, but we should use 
eqs. (15) and (16). 

Next, we compare our results with those by Kim and Kawasaki [28]. Their method is almost parallel to that we 
have used here. However, their basic equation is not FNH equations but Dean-Kawasaki equation. Thus, their MCT 
equation without interactions is the diffusion equation. On the other hand, our MCT equation (46) without memory 
kernel is equation for a damped oscillator. Thus, the existence of the momentum conservation equation in the basic 
equation naturally leads to the existence of the acceleration term in MCT equation. 

Third, let us compare our results with those by SDD [10], in which they used a simplified model of FNH. Although 
the approximation used here is similar to that used by SDD, the Galilean invariance is not preserved in their model 
equation. It implies that the violations of these conservation laws cause the artificial cutoff mechanism. 

Fourth, we compare ours with Das and Mazenko [8]. One of important differences between ours and theirs is that 
they regard V , which satisfies the constraint V = g/ p, as one collective variable. On the other hand, we introduce 
the new field variables 9 in eq. (15) and v in eq. (16) to satisfy FDR. Thus, their explicit expressions differ from ours. 
Second, let us discuss about their conclusion on the existence of the cutoff mechanism, i.e. Gpp oc Gpp = in the 
long time limit. As was indicated by ABL [23], the relation Gpp oc Gpp, used by Das and Mazenko (see eq. (6.62) 
in [8]), is not preserve FDR. Thus, we cannot conclude Gpp = from the relation Gpp = 0. However, the relation 
Gpp — might be valid, which was derived from their non-perturbative analysis. On the other hand, from eq. (A5), 
our FDR-preserving calculation under the first-loop order perturbation suggests Gpp = K{k)Gpp -\- Gep ^ 0, where we 
have used the numerical result of the non-ergodic parameter. Thus, our result in the first-loop order perturbation on 
Gpp is not consistent with Das and Mazenko. To resolve this discrepancy between their theory [8] and ours, or to verify 
their analysis on the cutoff mechanism, we need to find some identical relations without using the approximations. 
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B. Conclusion 



111 this paper, we reformulate a FDR-preserving field theory starting from FNH. When we assume that the correla- 
tions include the momentum decay fast enough, we have shown that there exists ENE transition and the noncrgodic 
parameter is equivalent to SMCT under the first-loop order approximation. These results give the theoretical basis 
of SMCT, and provides a route to go beyond SMCT. If we analyze correlations in higher-loop orders we believe that 
we will be able to construct a correct theory to explain the experimental and numerical results. 
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APPENDIX A: TIME REVERSAL TRANSFORMATION AND FDR 

From the linearity of the time reversal transformation, we can rewrite eq. (21) 
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When we express this time reversal transformation matrix as 0, eq. (21) can be represented by xp 
that the correlation functions satisfies 

G(-t) = OG{t) 0'^. 

Similarly, the self-energy satisfies 

where we have used the relation (30) or 
G-i = G„ 1 - Z. 



Oip. It implies 
(A2) 

(A3) 
(A4) 



Here, we avoid to write all the components of time reversal symmetry relations (A2) and (A3), because they are 
long and tedious equations. Instead, we present some typical relations, which are necessary for the calculation of 9p 
component of the SD equation; 



Gpp (fc, t) = Q i-t) [K (k) Gpp (fc, t) + Ggp (fc, t)) , 



(A5) 



G. {k,t)^-e{~t)dtGpp{k,t), 



G^^p (fc, t)^e i-t) (^G,„p (fc, t) + G,^p (fc, t) ) , 



G,^p{k,t)^-e{-t)dtGg^pik,t)^ 



(A6) 
(A7) 
(AS) 
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Here, we summarize some relevant relations among the self-energies 

s^^ (fc, t) = e (<) {dti:^^ (fe, t) - K (fc) E,-. (fc, t)) , (A9) 

E^,(fc,t) = -e(OE^.(fc,0, (Aio) 

s^^^ (fc, i) = -e (t) ( j^s^.^ (fc, t) - a^E,-,^ (fc, i)) , (All) 

5^e..(fc.*) = -0WS^§Jfe,i), (A12) 

E^^(fc,t) = See(fc>-0' (A13) 
^euAk,t) = -^^^^{k,-t). (AM) 

APPENDIX B: CALCULATION OF THE FDR-PRESERVING SD EQUATION 

In this APPENDIX, we calculate Op component of SD equation with the aid of APPENDIX A. From cqs. (23) 
and (31) , there is only —69 term in Sg which includes 6. Therefore, we obtain (Gg ^)g^(Ai — X2) = 5^b6{Xi — A2). 
Thus, 9p component of Go^^ ■ G satisfies 

[Go-i-G]^^(fc,t)=Gep(fc,i). (Bl) 
On the other hand. Op component of Z • G is expressed as 

[Z-G]g^(fe,t) = j ds\j:g^{k,t~s)Gpp{k,s) + J:gg{k,t~s)Gep{k,s) 

(fc, t - s) Gg^p (fc, s) + Eg-,^ (fc, t ~ s) G,^p (fc, s) 
+Eg-. (fc, t - s) Gpp (fc, s) + E^^ (fc, t - s) Gg^ (fc, s) 

+Eg.^ (fc, < - s) Gg„p (fc, s) + E^ (fc, t - s) Gs„p (fc, s) I . (B2) 
By using the relations (A9)-(A12), the first four terms in the right-hand side of eq. (B2) are rewritten as 
-J ds I (a.Egg (fc, t-s) + K (fc) Eg^ (fc, t - s)) Gpp (fc, s) + Eg^ (fc, t - s) Ggp (fc, s) 

+ [r^J^e-g^ (fe, ^ - + 5^5]^,^ (fc, t - s)^ Gg^p (fc, s) + E,-.^ (fc, t - s) G,^p (fc, s) 

~ (^7 0) Gpp (fc, i) 

- /" ds(E^^(fc,t-s)(if (fc)Gpp(fc,s) + Gep(fc,s))-Egg(fc,i-s)9,Gpp(fc,s) 
J —00 L 



^Eg.^ (fc,i-s) (^^Gg„p(fc,s) + G,„p(fc,s)^ -Eg^^ (fc,t-s)a,Gg„p(fc,s)|, 



(B3) 



where the boundary terms vanish except for T^gg {k,Q) Gpp {k,t) because of Eg. (fc,0) = from (A14) and 
Gpp (fc, -00) = Gg^p (fc, -00) = 0. 

Similarly, from the relations (A5)-(A8) the last four terms on right-hand side of eq. (B2) become 



/ ds \^gp{k, t - s) (K{k)Gpp{k, s) + Gep{k, s)) - Eg^(fc, t - s)dsGpp{k, s) 



^§a. (fe' i - (^^Gg^pik, s) + G,^p{k, s) ) - Eg. (fc, t - s) dsGg^p (fc, s) \. (B4) 
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From cqs. (B1)-(B4) , we obtain Op component of the FDR-preserving SD equation as 

Gep {k,t) + (fc, 0) Gpp ik,t) = - _^ { ^iJt^p * - (^) Gpp ^) + {k,s)) 

~'^gfj§ t — s) dgGpp (fe, s) — Sg^jp^ (fc, t — s) dsGg^p (fc, s) 

(fc,t-s) ('^Gg„p(fc,s)+G,„p(fc,s))|. (B5) 



Similarly, other components of SD equation can be obtained with the aid of the time reversal symmetry (Al). 

APPENDIX C: SOME EXACT RELATIONS OF EQUAL-TIME CORRELATION FUNCTIONS AND 

SELF-ENERGIES 

In this AppendiXjWe derive some relations for the equal-time correlation functions and the self-energies from the 
effective free energy F . Here, we note that the mean value is calculated by the canonical average over F. Since the 
equal-time correlation function satisfies 

= T {6p(r)K * Sp{r')) +T {6p{r)d{r')) . (CI) 
The Fourier transform of this equation becomes 

{Sp{k)e{^k)) = 0, (C2) 

where we have used the relation A'(fc) = Gpp{k,0). Thus, we obtain 

Gpe(fc,0)-0. (C3) 
Substituting (C3) into eq. (35), and setting i = 0, we obtain 

S^^(fc,0) = 0. (C4) 
Similarly, from the relations 

(^"(^^^^^'^J^) ^ TpoS^.Sir-r') 

= {go.ir)gp{r')) , (C5) 

and 

= {9.{r)9f{r')) +T {g^iry,{r')) , (C6) 

we obtain 

Gg„.,(fc,0) =G,,<,Jfc,0) = 0. (C7) 
Substituting eq. (C7), into eq. (37) at t = 0, we obtain 

Si>„p^(fc,0) = 0. (C8) 
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APPENDIX D: THE LIST OF THREE-POINT VERTICES 

In this Appendix, we write down the aU of the 4> included three-point vertices which are defined by (48) 



VppgAXi,X2,X:i) = -p-^\/r,J{X,-X2)S{X,-X3), (Dl) 

Vpp,AXl,X2,X3) = -TVr,a.S{Xi-X2)5iXi-X3), (D2) 

Vg (X,,X2,Xs) = ^—S{X,-X2)S{Xi-X3), (D3) 

V^^^^^{X,,X2,X3) - ~^J^t}HXi-X2)S{X^-X3), (D4) 

Vg^pp{Xi,X2,X3) = ~TS^d{Xi~X2)Vr,c.K{Xi-X3) 

+S{X,-X3)Vr,aK{Xi-X2)y (D5) 

Vg^pe{Xi,X2,X3) = -T6{Xi-X2)Vr,o.HXi-X3), (D6) 

^^3.9,9,(^1, ^2,^3) - -p^^l^5^pVrA5{Xi-X2)5{Xi-X:,)] 

+5p^5{Xi-X2)Vr^c.5{Xi-Xz)^, (D7) 

^ScS,3!^^ {.Xi,X2.X3) = -r|(5a/3Vri7[5 (^1 - X2) 5 {Xi - X3)] 

+ 6p^6{X^~X2)Vr^o.5{X^-X3)^, (D8) 

l/p„pg^ (Xi,X2,X3) = -^5o.p5{Xi-X2)5{X^-X:i), (D9) 



where (5(Xi - X2) = 5{ri - r2, h - ^2) and - X2) = 5{ti - t2)K * 5{rr - ra). 
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